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Sequential measurements on a single particle play an important role in experimental tests of the 
Kochen-Specker theorem and in Leggett-Garg inequalities. We provide a general method to analyze 
temporal quantum correlations, which allows to compute the maximal correlations in quantum 
mechanics. We present the full characterization of temporal correlations in the simplest Leggett- 
Garg scenario as well as in the most fundamental proof of the Kochen-Specker theorem. 
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Introduction — The physics of microscopic systems is 
governed by the laws of quantum mechanics and exhibits 
many features that are absent in the classical world. The 
best-known result showing such a difference is due to 
Bell pLJ. The assumptions of realism and locality lead 
to bounds on the correlations ~ the Bell inequalities - 
and these bounds are violated in quantum mechanics. 
Interestingly, this quantum violation is limited for many 
Bell inequalities and does not reach the maximal possi- 
ble value. For instance, the Bell inequality derived by 
Clauser, Horne, Shimony, and Holt (CHSH) bounds the 
correlations |2] 

{B) = {Ai®Bi) + {Ai®B2) + {A2®Bi)-{A2®B2), (1) 

where Ai and Bj are measurements on two different par- 
ticles. On the one hand, local realistic models obey the 
CHSH inequality {B) < 2, which is violated in quantum 
mechanics. On the other hand, the maximal quantum 
value is upper bounded by (B) < 2^2, a result known 
as Tsirelson's bound [3 . While this bound holds within 
quantum mechanics, it has turned out that hypothetical 
theories that reach the algebraic maximum (B) = 4 with- 
out allowing faster-than-light communication are possi- 
ble [Ij. This raises the question, whether the bounded 
quantum value can be derived on physical grounds from 
fundamental principles. Partial results are available and 
principles have been suggested that bound the correla- 
tions: In a world, where maximal correlations are ob- 
served, the communication complexity is trivial _5^, a 
principle established as information causality is violated 
[6], and there exists no reversible dynamics [7|. 

The question how and why quantum correlations are 
fundamentally limited has been discussed mainly in the 
scenario of bipartite and multipartite measurements. 
What happens, however, if we shift the attention from 
spatially separated measurements to temporally ordered 
measurements? There is no need to measure on distinct 
systems as in Eq. ([T]), but rather we may perform se- 
quential measurements on the same system. Then, an 
elementary property of quantum mechanics becomes im- 
portant: The measurement changes the state of the sys- 
tem. In fact, this allows to temporally "transmit" a cer- 
tain amount of information "5" and one would expect that 



the correlations in the temporal case can be larger than 
in the spatial situation. 

We stress that sequential measurements also have been 
considered in the analysis of the question how quan- 
tum and classical mechanics are different, the most well- 
established results here are quantum contextuality (the 
Kochen-Specker theorem [9]) and macrorealism (Leggett- 
Garg inequalities [lO]), cf. Fig. 1. The research in this 
fields has triggered experiments involving sequential mea- 
surements. For demonstrating such a contradiction be- 
tween classical and quantum physics, e.g. the correlation 

S5 ={AiA2)scq + (^2^3)scq + (^3^4)scq + (^4^5)scq 

- (A5Ai),eq (2) 

has been considered [TT1[T2]. Here, {AiAj)scq denotes a 
sequential expectation value, that is the average of the 
product of the value of the observables Ai and Aj when 
first Ai is measured, and afterwards Aj. One can show 
that for macrorealistic theories as well as for noncontex- 
tual models the bound {S^} < 3 holds, but in quantum 
mechanics, this can be violated. 

Here however, we are rather interested in the funda- 
mental bounds on the temporal quantum correlations. 
Special cases of this problem has been discussed before: 
For Leggett-Garg inequalities, maximal values for two- 
level systems have been derived [11] [13] and temporal 
inequalities similar to the CHSH inequality have been 
discussed [8^. 

We provide a method that allows to compute the max- 
imal achievable quantum value for an arbitrary inequal- 
ity and thus we solve the problem of bounding tempo- 
ral quantum correlations. First, we will discuss a simple 
method, which can be used for expressions as in Eq. ([2|, 
where only sequences of two measurements are consid- 
ered. Then, we introduce a general method which can 
be used for arbitrary sequential measurements, resulting 
in a complete characterization of the possible quantum 
values. Interestingly, our methods characterize temporal 
correlations exactly, while for the case of spatially sepa- 
rated measurements only converging approximations are 
known. 

Projective measurements — When determining the 
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FIG. 1: Sequential measurements occur in two different scenarios, (a) In the Leggett-Garg scenario one takes a single 
observable M which measures whether the physical system is in one of two possible macroscopic states. Then, one considers 
the correlations between these measurements at three different times, {M(ti)M(t2))seq, {M{ti)M{t3)}seq, and {M{t2)M{t3))scq- 
The values predicted by quantum mechanics contradict the assumption that the physical system is in any of these macroscopic 
states at any time and that the measurement reveals this state without disturbing it. (b) In the Kochen-Specker scenario one 
considers a set of observables Ai. Some of these observables are compatible and can therefore be measured simultaneously or 
in a sequence without any disturbance. Then one measures the correlations of simultaneous or sequential measurements of 
compatible observables, such as {Ai)scq, {AiAj)acq, and {AjAiAk)scq- For these correlations, one finds that quantum mechanics 
contradicts the assumption of noncontextuality. This assumption states that the result of a measurement should not depend on 
which other compatible observables are measured along with it. It should be noted, however, that the situation considered in 
this paper is more general than case (a) or (b), since no assumption about the time evolution or the compatibility of observables 
is made. 



maximal value for sequential measurements as in Eq. ^ 
we consider projective measurements, as these are the 
standard textbook examples of quantum measurements. 
The underlying formalism has been established by von 
Neumann |14l and Liiders [T5j. An observable A with 
possible results ±1 is described by two projectors n_|_ 
and n_ such that A — 11+ — n_ . If the observable A is 
measured, the quantum state is projected onto the space 
of the observed result, i.e. g i->- ^±gIl±/TT{gIl±). Apply- 
ing this scheme to the case of sequential measurements, 
one finds that the sequential mean value can be written 
as 

(A,A,),eq - ^[Tr{gAA,) + Tr{gA,A,)]. (3) 

It is interesting to notice that for pairs of ±l-valued ob- 
servables such a mean value does not depend on the order 
of the measurement [H]. 

The simplified method — We first show, how the max- 
imal quantum mechanical value for an expression like 55 
m Eq. ^ can be determined. First, we consider a set 
A = {Ai} of ±l-valued observables and a general ex- 
pression C = '^ij \j{AiAj)scq- The correlations given 
in Eq. ^ are just a special case of this scenario. Then, 
we consider the matrix built up by the sequential mean 
values Xij = {AiAj)seq- This matrix has the following 
properties: (i) it is real and symmetric, X — X"^ , (ii) 
the diagonal elements equal one, Xa = 1, and (iii) the 
matrix has no negative eigenvalue (or Xv > for any 
vector w), denoted as X > Q (see appendix). A similar 



construction for the matrix X, together with the opti- 
mization problem below, has been considered before in 
relation with Bell inequalities ;16J. However, our method 
involves a different notion of correlations, namely that 
given by Eq. (|3|. 

The main idea is now to optimize the expression 
C = '^ij over all matrices with the properties (i)- 

(iii) above. Hence, we consider the optimization problem 

maximize: Ay X^- (4) 

subjected to: X = X^ ^ and for all i, Xu = 1. 

Since all matrices X which can originate from a sequences 
of quantum measurements will be of this form, one per- 
forms the optimization over a potentially larger set. Thus 
the solution of this optimization is, in principle, just an 
upper bound on the maximal quantum value of 55. We 
will see that the bound is actually tight. Note that the 
optimization itself can be done efficiently and is assured 
to reach the global optimum since it represents a so-called 
semidefinite program |17j . In the case of this optimiza- 
tion can even be solved analytically and gives 

5,5 < ^(l + %/5) « 4.04. (5) 

It turns out that appropriately chosen measurements 
on a qubit already reach this value (see appendix and 
Refs. [HIIIH])- Hence, this upper bound is tight. More 
generally, one can prove that each matrix X with the 
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above properties has a sequential quantum representa- 
tion. Finally, note that if the observables in each se- 
quence are required to commute, then the maximal quan- 
tum value for 55 is known to be flqM = 4%/5 — 5 ~ 3.94 

The general method — The above method can only 
be used for correlations terms of sequences of at most 
two ±l-valued observables. In the following we discuss 
the conditions allowing a given probability distribution 
to be realized as sequences of measurements on a sin- 
gle quantum system in the general setting. We label as 
r = (ri, r2, . . . , r„) the results of an n- length sequence 
obtained by using the setting s = (si, S2, . . . , s„). The 
ordering is such that ri,si label the result and the set- 
tingfor the first measurement etc. The outcomes of any 
such sequence are sampled from the sequential condi- 
tional probability distribution 

P(r|s) = Pscq{ri,r2, . . • ,r„|si,S2, . . . ,Sn)- (6) 

In the case of projective quantum measurements each in- 
dividual result r of any setting s is associated with a 
projector 11*, which altogether satisfy two requirements: 
for each setting the operators must sum up to the iden- 
tity, i.e., X^r-n* = 1 and they satisfy the orthogonality 
relations II^II^, = Srr'^r^ where 5rr' is the Kronecker 
symbol. Finally, after the measurement with the setting 
s has given the result r, the quantum state is transformed 
according to the rule g U^QU^/P{r\s). 

In the following we say that the a conditional probabil- 
ity distribution P(r|s) has a sequential projective quan- 
tum representation if there exists a suitable set of such 
operators 11^ and an appropriate initial state g such that 

P(r|s) = Tr[n(r|s)n(r|s)tg], (7) 

with the shorthand n(r|s) = n^^n^^ . . .H'^^. 

Whether a given distribution P(r|s) indeed has such 
a representation can be answered via a so-called ma- 
trix of moments, which often appears in moment prob- 
lems [ini [^T^B5] . This matrix, denoted as M in the fol- 
lowing, contains the expectation value of the products of 
the above used operators n(r|s) at the respective posi- 
tion in the matrix. In order to identify this position wc 
use as a label the abstract operator sequence r|s for both 
row and column index. In this way the matrix is defined 
as 

Mr|,,r'|s' - (n(r|s)n(r'|s')^>- (8) 

Whenever this matrix is indeed given by a sequential pro- 
jective quantum representation, the matrix M satisfies 
two conditions: (a) Linear relations of the form Mi.|s;k|i — 
-^r'|s';k'|r if the underlying operators are equal as a con- 
sequence of the properties of normalization and orthog- 
onality of projectors, (6) M ^ since v^Mv > holds 
for any vector u, because such a product can be written 



(M(ti)M(t3)> 
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FIG. 2: Complete characterization of the possible quan- 
tum values for the simplest Leggett-Garg scenario. In this 
case, three different times are considered, resulting in three 
possible correlations (M(ti)M(t2)>Beq, {M{ti)M{t3))seq, and 
{M(t2)M(t3))seq. In this three-dimensional space, the pos- 
sible classical values form a tetrahedron, characterized by 
Eq. ||9| and variants thereof. The possible quantum mechan- 
ical values form a strictly larger set with curved boundaries. 

as the expectation value {CC^)g > which is non- 
negative for any operator C. Finally note that certain 
entries of this matrix are the given probability distribu- 
tion, for instance, at the diagonal Mr|s.i.|s = P(r|s). The 
main point, however, is the converse statement: Given 
a moment matrix with properties (a) and (h) above, the 
associated probability distribution P(r|s) always has a se- 
quential projective quantum representation. 

Hence, the search for quantum bounds represents again 
a semidefinite program. It is important to stress that this 
characterization is sufficient. This is in stark contrast 
with the analogue technique in the spatial Bell-type sce- 
nario [m 1^ , where one needs to use moment matrices of 
an increasing size n to generate better superset character- 
izations which only become sufficient in the limit n — > cx). 
However, indirectly, the sufficiency of our method has al- 
ready been proven in this context j22j (see appendix). 

Applications — In order to demonstrate the effective- 
ness of our approach, we discuss four examples. First, we 
consider the original Leggett-Garg inequality 

S = (M(ti)M(i2)>scq + (M(i2)Af(t3))«eq 

-(M(il)A/(i3)),eq< 1. ^' 

This bound holds for macrorealistic models and it has 
been shown that in quantum mechanics values up to 
5 = 3/2 can be observed [TUl El Our methods 

allow not only to prove that this value is optimal for any 
dimension and any measurement, we can for instance de- 
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termine all values in the three-dimensional space of tem- 
poral correlations {M(ti)M{tj)) which can originate from 
quantum mechanics. The detailed description is given in 
Fig. 2 and the calculations are given in the appendix. 

Second, we consider generalizations of the Eq. ^ with 
a larger number of measurements, known as iV-cycle in- 
equality [I1I2S], 



s 



N 



N-2 

E 

i=0 



2+1 /scq 



{An-iAq) 



(10) 



For this case, everything can be solved analytically (see 
appendix) leading to the bound 



Sn < N cos 



f- 



(11) 



which can be reached by suitably chosen measurements. 
This value has already occurred in the literature [TTJ [TH] , 
but only qubits have been considered. Our proof shows 
that it is valid in arbitrary dimension. Note that the fact 
that the maximal value is obtained on a qubit system is 
not trivial, although the measurements are dichotomic. 
For Kochen-Specker inequalites with dichotomic mea- 
surements examples are known, where the maximum 
value cannot be attained in a two-dimensional system 
[T5] and also for Bell inequalities this has been observed 

As a third application, we consider the noncontextu- 
ality scenario recently discovered by S. Yu and C. H. 
Oh |26j . There, thirteen measurements on a three- 
dimensional system are considered, and a noncontextu- 
ality inequality is constructed, which is violated by any 
quantum state. It has been shown that this scenario is 
the simplest situation where state-independent contextu- 
ality can be observed [57j , so it is of fundamental impor- 
tance. We can directly apply our method to the origi- 
nal inequality by Yu and Oh, as well as recent improve- 
ments [25] and compute the corresponding Tsirelson-like 
bounds. Our results are summarized in Table 1. 

Another class of inequalities is given by the guess- 
your-neighbor's-input (GYNI) inequalities [53], which, if 
viewed as multipartite inequalities, show no quantum vi- 
olation but a violation with the use of post-quantum no- 
signalling resources. We calculate the sequential bound 
for the case of measurement sequences of length three, 
instead of measurement on three parties. We consider 



P(000|000) + P(110|011) + P(011|101) 

+p(ioi|iio) < nc,Q <ns< n^s, 



(12) 



with the notation P(ri, r2, rajsi, S2, S3) as before, and 
possible results and settings G {0, 1} and Si G {0, 1}. 
We find that 



Ineq. 


NCHV 


State-independent 


Algebraic 


Sequential 




bound 


quantum value 


maximum 


bound 


Yu-Oh 


16 


50/3 ^ 16.67 


50 


17.794 


Opt2 


16 


52/3 ^ 17.33 


52 


20.287 


Opts 


25 


83/3 « 27.67 


65 


32.791 



1.0225, 



(13) 



TABLE I: Bounds on the quantum correlations for the 
Kochen-Specker inequalities in the most basic scenario. Three 
inequalities were investigated: First, the original inequal- 
ity proposed in Ref. [26) and the optimal inequalities from 
Ref. [28] with measurement sequences of length two (Opt2) 
and length three (Opt3). For each inequality, the follow- 
ing numbers are given: the maximum value for noncontex- 
tual hidden variable (NCHV) models, the state-independent 
quantum violation in three-dimensional systems, the algebraic 
maximum and the maximal value that can be attained in 
quantum mechanics for the sequential measurements. The 
latter bound is higher as the state independent quantum 
value, since the observables do not have to obey the com- 
patibility relations occurring in the Kochen-Specker theorem. 
Interestingly, in all cases the maximal quantum values are 
significantly below the algebraic maximum. 



while it is known that ^c,Q = 1 and ^Ins = 31 where 
the indices C, Q, S, NS label, respectively, the classi- 
cal, quantum, sequential and no-signalling bounds. So, 
in this case, the bound for sequential measurements 
is higher than the bound for spatially separated mea- 
surements. This also highlights the greater general- 
ity of our method in comparison with the results of 
Ref. [3130]: There, only temporal inequalities with se- 
quences of length two have been considered, where in 
addition the measurements can be split in two separate 
groups. In this case it turned out that the bounds were 
always reached with commuting observables. Our exam- 
ples show that this is usually not the case, when longer 
measurement sequences are considered. 

Discussion & Conclusions — For interpreting our re- 
sults, let us note that our scenario is more general than 
the scenarios considered by Leggett & Garg and Kochen 
& Specker. Leggett & Garg consider a special time evo- 
lution, g{t) — U {t) g{0)W (t) , which is mapped in the 
Heisenberg picture onto the observables. In our case, 
the observables are not connected via unitaries, this cor- 
responds to a more general time evolution. Compared 
with the Kochen-Specker scenario, our approach is more 
general since it does not assume that the measurements 
in a sequence are commuting. Nevertheless, if one wishes 
to connect existing noncontextuality inequalities to in- 
formation processing tasks, it is important to know the 
maximal quantum values (also if the observables do not 
commute), in order to characterize the largest quantum 
advantage possible. 

Furthermore, we emphasize that in our derivation it 
was assumed that the measurements are described by 
projective measurements and this condition is indeed im- 
portant. In fact, this sheds light on the role of projective 
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measurements: One can easily construct classical devices 
with a memory, which give for sequential measurements 
as in Eq. ([2| the algebraic maximum S5 = 5. These clas- 
sical devices must also have a quantum mechanical de- 
scription. Our results show, however, that in this quan- 
tum mechanical description more general than projective 
measurements must occur and a more general dynamical 
evolution than the projection is required. From this per- 
spective, our results prove that the memory that can be 
encrypted in quantum systems by projective measure- 
ments is bounded. 

Our results lead to the question, why quantum me- 
chanics does not allow to reach the algebraic maximum 
of temporal correlations, as long as projective measure- 
ments are considered. We believe that proper generaliza- 
tions of concepts like information causality and commu- 
nication complexity might play a role here, but we leave 
this question for further research. A first step in explain- 
ing quantum mechanics from information theoretic prin- 
ciples lies in the precise characterization of all possible 
temporal quantum correlations, and our work presents 
an operational solution to this problem. 

Acknowledgements — We thank J.-D. Bancal, T. Fritz, 
Y.-C. Liang, and G. Morchio for discussions. This 
work has been supported by the EU (Marie Curie GIG 
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Appendix 

Al: Discussion of the simplest Leggett-Garg 
scenario 

In this part we provide some further details how to 
determine the set of possible quantum values for the 
simplest non-trivial Leggett-Garg scenario as shown in 
Fig. [2j Here it is assumed that one can measure an ob- 
servable M at three different time instances ^1,^2 7^3 as 
shown in Fig. [T] which gives rise to three different ob- 
servables Ai = M{ti) with i = 1, 2, 3. 

However, rather than being interested in determin- 
ing the full sequential probability /'(rjs) for all possible 
combinations we are here only interested in some lim- 
ited information, namely only for the correlation space. 
This means that from a general distribution we only 
want to reproduce the correlations terms {AiAj)sc(i with 
1 < * < j < 3 each defined by 

(A,A,),eq = P{n = r,\i,j) ~ P(r, ^ r,\i,j). (14) 

Thus we want to characterize the set 

^qm — {<lij G R"^ '■ Qij = (^z^j)seqj 

{AiAj)scq has projective quantum rep.}. (15) 



For this we refer to problem given by Eq. Q, with 



1 (^1^2)scq (^1^3) scq 

scq 

(^1^3)scq (^2^3)scq 1 



(16) 



Now we know that any matrix of this form has a sequen- 
tial projective quantum representation if and only if X is 
positive semidefinite. However a matrix satisfies A )^ 
if and only if the determinant of all principles minors are 
non-negative. This finally gives 



^qm — {qij € IR'^ 



< 1, 



1 -I- 2(712913923 > 9l2 



2 

lis 



923}- 



(17) 



which is the plotted region of Fig. [2] 

At last let us comment that via the general method 
one can also determine the probability distribution of 
a general scenario. However, since the general method 
requires the solution of a SDP, where in particular some 
entries are completely unknown, it is not so easy to write 
down a simple condition on the correlations. 



A2: Detailed discussion of bounds for the A-cycle 
inequalities 



We first need the general form 20J for Eq. ( 10 1 



(18) 



i=0 



with 7 = (70, . . . , 7Ar-i) G {^li 1} and an odd number 



of —1. Every expression in Eq. (18 1 can be obtained 
from another one [in particular, from Eq. ( 10 )] by the 
substitution Ai — > —Ai for some indices i. Since the 
quantum bound does not depend on the particular 7, for 
the case odd N, we can consider the expression 



N-l 



seqj 



(19) 



with index i taken modulo N. The optimization problem 
in Eq. Q, therefore, can be expressed as 



maximize: -Tr{WX) 
subjected to: X >z and Xa — 1 for all i 
where W is the circulant symmetric matrix 

1 ... 1 

10 1 



(20) 



W 



1 





1 



■• 1 

1 



(21) 
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The condition X ^ 0, i.e. Xv > for any real vector and the maximization problem can be expressed as a SDP 
V, follows from the fact that {AiAj)scq = 5Tr[gi(AiAj + 
AjAj)] and the fact that the matrix Y — Tr[g{AiAj)] 



fulfils V Yv > for any real vector v, and X is the real 
part of Y. 

By using the vector A ~ (Ai, . . . ,Xn), the dual problem 



for the semidefinite program in Eq. ( 20 ) can be written 



as (see Ref. |J/7j for a general treatment and Ref. [TB] for 
the discussion of a similar problem) 



minimize: TT{diag{X)) 
1 



subjected to: 



W + diag{X) h 



(22) 



where diag{\) denotes the diagonal matrix with entries 
Ai, . . . , Aat. 

Let us denote with p and d optimal values for, respec- 
tively, the primal and the dual problem, then d > p. We 
shall provide a feasible solution for the dual problem with 
d = N cos(^) and a feasible solution for the primal prob- 
lem with p ^ d, this will guarantee the optimality of our 
primal solution. 

We start by finding the maximum eigenvalue for W. 
Since is a circulant matrix, its eigenvalues can be writ- 
ten as ED 



-2cos,'^) 



(23) 



1, and the maximum is 



for j = 0,. 

f^max — 2 cos ( ) • 

For a pair of Hermitian matrices A, B, it holds 

fJ-miniA + B) > IJrmniA) + llynin{B) , wherC p,rain 

denotes the minimum eigenvalue. Therefore, A = 
(cos (-^) , . . . , cos (-^)) is a feasible solution for the 
dual problem and Tr[diag(A)] = A^cos(-^), and p < 
iVcos(^). 

Now consider the matrix X'^^j — {xi,Xj), with 
, a; TV unit vectors in a 2-dimensional space such 



Xi, 



that the angle between Xi and Xi^i is ^^-^t^, and (•, •) de- 



noting the scalar product. Clearly, X' is positive semidef- 

(^)forz = 0,...,iV-l 



inite, X-^ = 1 and X'^ = — cos i 
with sum modulo N . 

It follows that p ^ d = N cos (-^) and the solution X' 
is optimal. 

To prove that X' can be obtained as matrix of expec- 
tation values for sequential measurements, it is sufficient 
to notice that given two 3-dimensional unit vectors a, 6, 
and defined <Ja = a ■ a and at = a • b, then, by Eq. (Isj), 
{o'aO'b)seq = a ■ b, independently of the initial quantum 
state Q. In fact, explicit observables reaching this bound 
have already been discussed in the literature [TTJ HH] . 

For the case N even, we can consider the expression 



N-2 

E 

i=0 



{AiAi 



{AqAn^i), 



(24) 



as in Eq. (20 1, with the proper choice of the matrix W. 
Such a SDP has been solved in Ref. [16]. The solution is 
analogous to the previous one: A set of observables, for 
a two-level system, saturating the bound, again, inde- 
pendently of the quantum state, is given by observables 
Ai — d ■ Xi, where the vectors Xi are on a plane with an 
angle separating Xi and Xi+i. 

As opposed to the N odd case, such a bound can be 
also reached with commuting operators, this corresponds 
to the well known maximal violation of Braunstein- Caves 
inequalities |T^. 

The above results prove that the bound computed in 
Ref. |18j for sequential measurements on qubits, coin- 
ciding with the value explicitly obtained in Ref. [TT], is 
valid for any dimension of the quantum system on which 
measurements are performed. 

Finally, we stress that the construction of the above 
set of observables from the solution of the SDP, i.e. the 
matrix X or the set of vectors {xi} such that Xij = 
{xi,Xj), is actually completely general. We recall that 
the vectors {xi} can be obtained, e.g., as the columns 
of the matrix X^^^, which is well defined since X 0, 
and, therefore, the dimension of the subspace spanned by 
them is equal to the rank of the matrix X. In the previous 
case, since we were dealing with vectors in dimension 
d < 3, we used the property of Pauli matrices 



{o-ajCTfc} = (TaCTb + Cr^aa = 2(5- 6)1. 



(25) 



For matrices X with higher rank, the corresponding 
vectors {xi} will span a real vector space V of dimension 
d > 3. Now for general complex vector spaces V with a 
symmetric bilinear form ( , ), an analogous of Eq.(25l, 
namely 



{A^,, A„} = 2(w,it)l, foranyw,t;ey (26) 

defines the associated Clifford algebra, where Au , Ay are 
the associated elements of the algebra and 1 is the iden- 
tity. In particular, for V a real vector spaces with 
the standard Euclidean scalar product, such a prop- 
erty is satisfied by the corresponding canonical anti- 
commutation relations (CAR) algebra. The latter has 
always a representation in terms of self-adjoint opera- 
tors on an Hilbert space, e.g., operators written in terms 
of Fermionic creation and annihilation operators on an 
antisymmetric Fock space |32j . or written in terms of 
the Clifford algebra generators explicitly constructed in 
Tsirelson's original argument [33l [34] 

As a consequence, for every positive semidefinite real 
matrix X with diagonal elements equal to 1 , one can find 
a set of unit vectors {xi} giving Xij = (xi, Xj) and a set of 
±l-valued observables {Ai}, associated with {xi} , such 
that 



{AiAj) seq — Tr 



giA.Aj + AjA,) 



(x„x,), (27) 
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for all quantum states q. In particular, if the rank of X is 
d, such operators can be chosen as 2^* x 2"* Hermitian ma- 
trices [32] . This shows the completeness of the simplified 
method. 



A3: Completeness of the general method 

In this part we shortly comment on the completeness 
of the presented general method. As pointed out, this 
has already been proven indirectly in the context of the 
spatial bipartite case [22]. 

At first let us change slightly the notation in order to 
make it closer to the one used in Ref . [22] . In the follow- 
ing we do not explicitly consider the matrix M from the 
main text, but rather a slightly smaller matrix where one 
erases some trivial constraints. In the following the set 
{Ei} contains all projectors 11^, but one of the outcomes 
k from each setting s is left out. Note that we also use 
a single subscript to identify setting and outcome. Then 
one looks at the matrix 



Xl^ = TT[E{n)E{vy p] 



(28) 



being built up by all products u — (ui, U2, . . . , w;), 
E{vl) = EuiEu2 ■ ■ ■ Eu, of the operators {Ei} of at most 
length / < n, and the single extra "sequence" u = of 
the identity operator E{0) = 1. Again this matrix has 
to satisfy hnear relations parsed as xZv — Xu'v'; 
operators fulfill E{u)E{vy — E{u')E{'v'y as a conse- 
quence of the orthogonality properties of projectors, and 
that x" h 0. 

That this matrix is positive semidefinite can be ver- 
ified as follows: Let us first assume that there exists a 
sequential projective quantum representation. Consider 
the operator C = J^u Cu-E'(u)^' with arbitrary Cu € C! and 
evaluate the expectation value of CC^, which provides 

Ti{CC^g) = ^ CuTr[E{u)'fE{v)g]cl (29) 

U.V 

= Ec"^uv<>0 (30) 

u,v 

using linearity of the trace. The final inequality holds be- 
cause CC^ >: and p ^ are both positive semidefinite 
operators. Since Cu G (D are arbitrary the condition given 



by Eq. (30) means that x" ^ is positive semidefinite. 



For the reverse one needs a way to construct an ex- 
plicit sequential projective quantum representation out of 
the matrix xCv satisfying the above properties. For this, 
clearly more difficult part, we refer to Ref. [52] and just 
mention the solution. For the given positive semidefinite 
matrix x" one associates a set of vectors {|eu)} by the 
relation Xuv = (eu|ev)- From this set of vectors one now 
constructs an appropriate state and corresponding pro- 
jective measurements hy H = span({|eu)}), p = |eo)(eo|, 
and Ei = proj(span({|eu) : ui = i})) where proj means 



the projector onto the given subspace. That these solu- 
tion satisfies all the required constraints are shown in the 
proof of Theorem 8 of Ref. [22]. 

However, in the spatial case considered in Ref. [22], 
some of these operators, additionally, have to commute 
since they should correspond to measurement onto dif- 
ferent local parts. This cannot be inferred, in general, 
by a finite level x" and this is the true reason why their 
characterization gets only exact (for the commutative bi- 
partite Tsirelson problem) if n tends to infinity. How- 
ever, since we do not require here that the measurements 
of different settings have to commute we, luckily, do not 
need this. 
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